Lifetimes of atoms trapped in an optical lattice in proximity of a surface 
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We study the lifetime of an atom trapped in an optical vertical lattice in proximity of a massive 
surface using a complex scaling approach. We analyze how the presence of the surface modifies the 
known lifetimes of Wannier-Stark states associated to Landau-Zener tunnelling. We also investigate 
how the existence of a hypothetical short-distance deviation from Newton's gravitational law could 
affect these lifetimes. Our study is relevant in order to discuss the feasibility of any atomic- inter- 
ferometry experiment performed near a surface. Finally, the difficulties encountered in applying the 
complex-scaling approach to the atom-surface Casimir-Polder interaction are addressed. 
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I. INTRODUCTION 

The problem of atoms trapped in an accelerated peri- 
odic lattice is an old one in solid state physics [TH51 . How- 
ever, some phenomena like pseudo-eigenstates (Wannier- 
Stark ladders) and the lifetime limitations of those states 
(Landau-Zener tunneling) are still a current subject of 
research both on the experimental side and theoretically 
[51[T5] (and see [5] and reference therein for a more exten- 
sive review). In this work, we study the Landau-Zener 
effect [TH - fTB] for atoms trapped in a vertical standing 
wave in the presence of short scale interactions. 

This is done in the context of the theoretical mod- 
eling of the experiment FORCA-G (FORce de CAsimir 
et Gravitation a courte distance) aiming at measuring 
the Casimir-Polder force acting on a neutral atom in the 
presence of a massive dielectric surface by atomic interfer- 
ometry. At the same time, the experiment is configured 
to search for non-Newtonian deviations from the gravi- 
tational law at short distances [T7l - f23] . These two goals 
are reached by the use of atoms trapped in a vertical 
optical standing wave in front of a massive surface real- 
izing the trap. Such a lattice allows a precise knowledge 
of the distance between the atom and the surface, this 
distance being a multiple of the wavelength of the trap. 
The atomic states in such a trap were the main subject 
of a previous paper [19] . In that paper the shape of the 
wavefunctions and the modification of energy levels due 
to the presence of the surface were deduced by solving a 
standard eigenvalue problem. In this context, the states 
are supposed to be bound states of the system, meaning 
that their lifetime is infinite. 

Nevertheless, solid state physics teaches us that the 
states in an accelerated periodic lattice are not bound 
states but resonance states with a finite lifetime. The 
subject of this paper is to quantify these lifetimes for 
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the states calculated in [19j in order to check the feasi- 
bility of the experiment, knowing that the measurement 
of the interaction by atom interferometry takes at most 
a few seconds. We want to understand, more in detail, 
how the presence of the surface and of a hypothetical 
non-Newtonian gravitational potential modifies the well- 
known finite lifetimes of standard Wannier-Stark states. 
This paper is organized as follows. In Sec. |TT] we de- 
scribe the system under scrutiny and the method we use 
to calculate the lifetime of the atomic states. Section ITlll 
is devoted to the effect of the surface. This section is sep- 
arated in two parts. As a first step, we study the effect 
of a surface treated as a boundary condition for an atom 
below it. In the following, we investigate the phenomena 
arising when tunneling of the atom through the surface 
is allowed. This is done for an atom above the surface 
in order to clearly discriminate the effect of the finite- 
ness of the potential barrier from other resonance effects. 
Then, Sec. |IV| shows the effect of a deviation from New- 
ton's gravitational law at short distance on the lifetimes 
of the trapped atom whereas Sec. |V] discusses the prob- 
lem arising when the Casimir-Polder interaction is taken 
into account. In this part, the difficulties of the treat- 
ment of the Casimir-Polder effect by the complex-scaling 
method are highlighted. 



II. LIFETIMES OF ATOMS IN 
WANNIER-STARK STATES 

Let us consider the Hamiltonian of an atom of mass rUa 
trapped in a periodic lattice V{z) = ^ (1 — cos{2kiz)) in 
the presence of a linear gravitational potential 
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where U is the depth of the trapping potential, ki repre- 
sents the wavevector of the trap and g is the earth gravi- 
tationnal constant. This operator is the sum of the well- 



known Block Hamiltonian Hg 
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[T] and a linear potential —niagz (the z axis is oriented 
downwards). As a consequence, the eigenstates of Eq. 
(fTl) are derived from the Bloch states, generally taking 
into account one single Bloch band. The spectrum ob- 
tained using this procedure can be written as 
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where ^ is the periodicity of the lattice, A; being its 
wavelength, Sa is the mean value of the energy of the 
Bloch band a under consideration and n is an integer 
labelling the well where the state is centered. The states 
obtained are known as Wannier-Stark states and were 
first discussed for electrons in a crystal submitted to a 
constant electric field ,2J. The existence of such eigen- 
states results from the single-band approximation. Nev- 
ertheless, a calculation ignoring this approximation con- 
firms the existence of these states, but associates a finite 
lifetime to each of them. This lifetime physically corre- 
sponds to the so-called Landau-Zener tunneling, namely 
the possibility for a particle in a Wannier-Stark lattice 
site to transist from one Bloch band to another TltilGj . 
In other words, the discrete spectrum in Eq. (pi) is im- 
mersed in the continuum of eigenvalues of the Hamilto- 
nian (IT) and can thus be seen as a resonance spectrum 
[Mj . The spectrum of Eq. ([2| should thus be rewritten 
as [3H5] 
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The eigenstates associated to the spectrum (Is]) are 
metastable states with a finite lifetime r = jl^. In or- 
der to completely characterize the states of an atom in 
our optical lattice, we need to work out this complex 
spectrum to evaluate their lifetime in the trap. 

Several methods have been developed so far in order 
to study this complex spectrum pSHSOj . each using ap- 
proximations such as a periodic potential with a finite 
number of gaps or a finite periodic lattice. Most of these 
methods are based on the analogy of the Wannier-Stark 
sates with scattering states. In this theoretical frame- 
work, the complex spectrum (^1) can be calculated as the 
set of poles of the scattering matrix of the system. The 
method we will use in this work is based on the rota- 
tion of the Hamiltonian M in the complex plane and the 
derivation of the complex spectrum as eigenvalues of the 
obtained non-hermitian Hamiltonian. This approach is 
known as complex- scaling method [31 1311 132] and the ro- 
tation of the Hamiltonian is performed via the operator 
C{e) such that 



C{0)ip{z) = <p{ze" 
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The one-dimensional complex-scaled Hamiltonian is then 
H^l = "-^ + I (1 - cos{2kize^')) - magze'. (5) 
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FIG. 1: Schematic representation of the complex eigenvalues 
of the complex-coordinate scaled Hamiltonian H^^' according 
to the Balslev-Combes theorem 1331. 



According to the Balslev-Combes theorem [351 , the com- 
plex spectrum ([3| can then be found by solving the eigen- 
value problem 



H^^f^^(t>a,n) = E, 
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where the real part of a given eigenvalue represents its 
energy level, while its imaginary part is the width of the 
resonance. 

We stress here that resonance wavefunctions are char- 
acterized by a divergent behavior when z — >■ cx) as 
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k being the wavevector of the atom and S{k) the scat- 
tering matrix of the system. As demonstrated in [31 , the 
divergent part of the resonant wavefunctions in Eq. ([TJ) 
is regularized by the complex scaling. This is not the 
case for a non-resonant wavefunction of the continuum. 
As a consequence, it is important to distinguish the reso- 
nant from the non-resonant part of the complex spectrum 
worked out solving Eq. (rol. A schematic representation 
of the eigenvalues of Hamiltonian ^ is shown in Fig. [l] 
We see in this figure that the continuum eigenvalues lie on 
a line forming an angle 26 with the real axis, being the 
angle chosen for the transformation ^ . At the same time 
the real bound states are located on the real axis (corre- 
sponding to the fact that they have an infinite lifetime) 
while the resonance states are the complex eigenvalues 
which are not part of the rotated continuum. In the case 
of the Wannier-Stark Hamiltonian, the presence of en- 
ergy bands due to the periodicity of the trap leads to a 
peculiar structure of the complex eigenvalues. The eigen- 
values of the Hamiltonian (Isl) are shown in Fig. [2] These 
eigenvalues are obtained using the pseudo-periodicity of 
the trap. Indeed, the fact that we calculate states in an 
infinite pseudo-periodic lattice allows us to restrict the 
search of the spectrum in Eq. (Isl to the analysis of the 
spectral properties of the Floquet-Bloch operator [S] 
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over a Bloch period [H [S]. The eigenvalues are then 
worked out using a Discrete Variable Representation 
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FIG. 2: Complex eigenvalues for the 3 first Bloch bands of 

(9) 

the Hamiltonian H^g. The energies are expressed in units 



trap with wavevector ki (in the case of FORCA-G, Er — 5.37- 
10"^" J and r = lEr corresponds to r = h/T = 1.96 • 10"^ s). 
The spectrum here is calculated for u = 3Er. 



(DVR) method and a QR algorithm for complex sym- 
metric matrices [M] . 

In order to check the validity of our method, we com- 
pare the real parts of the calculated eigenvalues with 
the spectrum obtained using two different methods. The 
first is the tight-binding model O [17] , a standard semi- 
analytical resolution of the Schrodinger equation for -ffws 
which does not give access to the finite lifetime of the 
states. The second one is a finite-difference numerical ap- 
proach to -ffws without complex scaling. Table |l] shows 
the good agreement between the real parts of eigenvalues 
^ and the second approach. The difference between our 
results and the tight-binding model can be traced back 
to the fact that this approach restricts the Hilbert space 
to the first Bloch band. 



Concerning the imaginary parts, representing the life- 
time of the band under scrutiny, Niu et al. have proposed 
a general formula to evaluate the inverse lifetime of the 
standard Wannier-Stark states [321 ■ This formula gives 
an estimate of the tunneling rate 7 as a function of a crit- 
ical acceleration ac and the dimensionless accelaration a 
due to the linear potential. 
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In our case, a — Sri and the critical acceleration can be 



expressed as 
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where A is the half width of the energy gap between the 
first and the second Bloch band and _ftr = ^ is the wave 
number of Bragg scattering corresponding to the n-th 
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1.22663 
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1.15645 


0.15656 
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1.08638 


1.08638 


0.08649 



TABLE I: Table of the real parts Re(i5i,„) of the complex 
energies ([3| compared with the Wannier-Stark spectrum £i^„ 
obtained numerically without complex scaling and with the 
energies £n resulting from a semi-analytical calculation in the 

first Bloch band. The energies are given in units of Er = -2^ 
which is the natural unit for energy levels in a trap with a 
wavevector ki. The numerical example corresponds to f/ = 

SEr. 



gap (here we will take n — 1). With this expression, we 
obtain, for the second Bloch band, a tunneling rate of 
the order of 0.003 i^^ which is of the same order of mag- 
nitude of the value of the complex part of the spectrum 
([3]) as it can be seen on Fig. [2] An additional verifica- 
tion with other values of the well depth assures us that 
the complex-scaling method provides satisfying results 
for this physical system. As a consequence, in the next 
section, we will use this method to analyze the lifetimes 
of the Wannier-Stark states in the presence of a surface, 
whose main features have been already discussed in [19| . 



III. LIFETIMES IN THE PRESENCE OF THE 
SURFACE 

A. Presence of the surface 

In our problem, we consider the presence of a surface 
at z = 0. Our atom is thus located below the surface, 
in the region 2: > (we recall that the z axis is oriented 
downwards). The presence of this boundary condition 
breaks the quasi-periodicity of the system. The poten- 
tial modifying the optical trap is no longer linear, since 
it must be considered as the gravitational linear poten- 
tial for z > and an infinite potential barrier for z < 0, 
describing the impossibility of the particle to penetrate 
the mirror. We have shown in a previous paper jl9| that 
this surface induces a modification of the energies and 
states of the Hamiltonian i?ws- We now want to ver- 
ify that this modification does not reduce drastically the 
lifetimes of our atoms in the trap. 

As anticipated before, the presence of a boundary con- 
dition in z = breaks the translational symmetry of 
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FIG. 3: (Color online) Lifetimes of the first Bloch band in 
function of the quantum number n labelling the well of the 
trap for five different depths. From top to bottom we have: 
U = 2Er (blue dots), U = l.^Er (green dots), U = lEr 
(orange dots), U = Q.bEr (red dots), U — Q.lEr (black dots). 
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FIG. 4: (Color online) Lifetime of the first Bloch band calcu- 
lated with the complex scaling method (Red dots) in function 
of the trap depth. The calculated lifetimes are compared with 
the Landau-Zener formula (Black line) for standard Wannier- 
Stark sates. 



the problem. Thus, we can no longer use the quasi- 
periodic approximation and we have to solve directly 
the eigenvalues problem ([6]). This is done using a finite- 
element method and once again a QR algorithm for sym- 
metric complex matrices using the subroutine zgeev of 
the numerical package LAPACK 36 . However, using this 
method, numerical problems arise when the well depth 
is larger than U = 2.2Er, since in this case the imagi- 
nary parts of the complex eigenvalues are too small. As 
a matter of fact, the lifetime of a given Wannier-Stark 
well is a strongly increasing function of the well depth. 
This dependence can be well understood using the image 
of tunneling. Indeed, the finite lifetime of the metastable 
states can be seen as a consequence of the possibility for 
an atom in a lattice site n to tunnel through the poten- 
tial toward an upper Bloch band in another site. This 
is favorized by the linear potential which induces reso- 
nant tunneling between the first band in a given site to 
e.g. the second band in a farther site. As a consequence 
the atom reaches a new state with a new energy destroy- 
ing the former one. Moreover, the atom can reach a state 
with an energy higher than the trap depth, being thus no 
longer trapped. So, the tunneling probability depends on 
the lattice depth: the deeper the trap is, the weaker this 
probability becomes. 

This feature is illustrated in Fig. [3) where lifetimes 
are calculated for different well depths. Here we have 
taken into account only the lifetime of the first Bloch 
band. Indeed, as it can be seen in Fig. [2| the lifetime 
of the first band is the longest one, and the larger the 
band under scrutiny is, the shorter is this lifetime. In 
the particular case of our experimental parameters, the 
second Bloch band has a mean energy of 62 ~ 5.45 Er 
and is then not trapped by a well depth U = SEr. As 
a consequence, only the lifetime of the first band is of 
interest for us. 



The analysis of Fig. [3] shows, moreover, that for any 
well depth, the first-band lifetimes are function of the 
considered well. In particular, the lifetimes do not vary 
remarkably starting approximately 5 wells away from the 
surface, whilst being functions of the considered well in 
proximity of the boundary z = 0. In order to verify 
the coherence of the calculation the constant lifetimes far 
from the surface (function of the well depth) can be com- 
pared to the ones obtained through the Landau-Zener 
formula, derived for the standard Wannier-Stark states 
(i.e. in an infinite lattice). This comparison is shown in 
Fig. [4] In this figure, we can see that the Landau-Zener 
formula still gives a good description far from the surface. 
This suggests that the lifetimes of our metastable states 
in the presence of a surface are of the order of 10^'* s for a 
well depth oiU — 3Er. However, in this analysis we have 
treated the mirror as an infinite barrier at z = 0, which is 
not very realistic. It is of interest to investigate how the 
results are modified if the mirror is modeled by means 
of a finite potential barrier. This will be done assuming 
that the atom is located above the mirror. 



B. Atom above the mirror 

A way to increase the lifetime of the states is to place 
the atom above the surface rather than below. In this 
case, if we consider as before the mirror as an infinite 
barrier, the potential is a well of infinite height. As a 
consequence, it supports only a discrete set of real eigen- 
states having an infinite lifetime. In order to model the 
surface more precisely, we allow the tunneling of the atom 
into this surface. Therefore we replace the infinite bar- 
rier by a finite one as shown in Fig. [51 As a consequence, 
the complex-scaled Hamiltonian of the system can be ex- 




where Vq € R. In order to study the behavior of the 
atoms in front of a surface with a finite size, we in- 
vestigate the complex spectrum for a given well depth 
([/ = lEr in the numerical example) as a function of 
the height of the barrier chosen to represent the surface. 
More specifically, we focus our attention on the transi- 
tion from bound states to resonances for atoms in front 
of a potential barrier with a varying height. 
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FIG. 5: Shape of the potential for atoms above a finite surface. 
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The resulting complex eigenvalues for different heights 
of the barrier are shown in Fig. [6J 
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FIG. 6: (Color online) Complex spectrum of the Hamiltonian ( |11[ | describing an atom in front of a potential barrier with 
different heights (from left to right and top to bottom Vb ~ Cl-Br, Vb = lEr, Vo = 2Er and Vb — 5Er) (black dots) in 
comparison with the standard Wannier-Stark lifetimes (blue lower horizontal line). The depth of the periodic trap is chosen 
to he U = Er- The solid red fine represent the line of the continuum with a slope of 28 and the horizontal green upper line 
highlights the eigenvalues with an imaginary part equals to 0. 



In this figure, we see that the presence of the barrier induces the emergence of a continuum associated with 



the standard resonances spectrum. The position of this 
continuum depends on the height of the barrier: more 
precisely, the continuum region of the spectrum always 
starts at a value around Vq. When the barrier is very 
low (i.e. of the order of the depth of the periodic trap 
or below) , we observe that the lifetime of the atom close 
of this barrier is shorter than the lifetime at longer dis- 
tances. This can be understood thanks to the interpreta- 
tion of Landau and Zener of the resonance phenomenon. 
Indeed, if we consider that the finite lifetime of the atom 
arises from the resonant tunnelling through the trap, it is 
clear that this tunnelling has a higher probability when 
the atom is close to the barrier because it is not trapped 
on the left side. In addition, we observe that when the 
barrier becomes higher than the well depth, some bound 
states appear corresponding to the closest states submit- 
ted to a potential well. 

As expected, the effect of a finite barrier representing 
the surface is completely negligible for the wells far from 
the surface and starts to be visible for the closest wells 
when the barrier height is of the order of the well depth or 
below. This means that the penetrability of the surface 
should be comparable to that of the periodic trap in order 
to allow a significant modification of the behavior of the 
complex spectrum. This is reasonably not at all the case 
for a solid mirror and we can conclude that the finiteness 
of the surface should not play a major role in the lifetime 
of our states in the trap. 



IV. NON-NEWTONIAN GRAVITATION 

After observing that the presence of the surface 
does not reduce drastically the lifetime of our modi- 
fied Wannier-Stark states, we have to verify that this 
is not the case even in the presence of a Yukawa gravita- 
tional potential. Indeed, the main goal of the experiment 
FORCA-G is to search for a hypothetical deviation from 
Newton's law at short distance predicted by some unifica- 
tion theories. We have then to estimate the modification 
of the lifetimes of the states in the presence of such a de- 
viation which can be written as an additional potential 
of the form [H] 



Uy = 2TTaYGpsmaXYe 



-2x/Ay 



(12) 



where G = 6.67 x 10" 



'kg- 



is the universal 



gravitational constant and ps is the density of our sur- 
face (for which we have chosen the density of the sili- 
con Ps = 2.33 X lO'^kgm"^). We want to stress here 
that the Yukawa potential presented in Eq. (12 1 de- 
pends on two parameters ay and Ay- These parameters 
represent respectively the coupling strength of the non- 
Newtonian deviation and its typical range. The aim of 
experiments devoted to non-Newtonian gravitation is to 
impose constraints on the value of these two parameters. 
The present constraints and the ones predicted for the 
experiment FORCA-G are shown in Fig. [7] 
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FIG. 7; (Color online) Present constraints (yellow region) and 
predicted constraints (coloured lines) for the FORCA-G ex- 
periment for the Yukawa potential in the (ay, Ay) plane. The 
colored zones are excluded. This figure is taken from |19| and 
adapted from [37] . 



In order to analyze the modification to the lifetimes 
induced by a hypothetical Yukawa term, we have to cal- 
culate the spectrum of the Hamiltonian 

H —Hws + Uy 

p^ U 

= - h — (1 - cos(2fc;z)) - mapz (13) 

zrua 2 

+ 2TTayGpsmaXye-^'-/^^ . 

As done previously, we apply the transformation Q to 
this Hamiltonian, obtaining the non-hcrmitian Hamilto- 
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The complex spectrum resulting from this modified 
Hamiltonian is presented in Fig. IS] This figure shows 
that even if the presence of a deviation from Newton's 
gravitation law modifies the lifetimes of the atom in the 
trap, particularly for the closest wells, this lifetime re- 
mains very large in comparison with the duration of the 
experiment. Indeed, for practical reasons due to the co- 
herence of the cold atoms in the experimental setup, the 
duration of the measurement in FORCA-G will not ex- 
ceed Is. As a consequence, the Landau-Zener effect is 
not a limitation for the measurement of a deviation from 
Newton's law in our experiment. This can be more ac- 
curately verified from the analysis of Table III) 



This table confirms that the lifetimes of an atom in the 
trap are much longer than the duration of the measure- 
ment estimated in |22j. 
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FIG. 8; (Color online) Values of the complex spectrum in the 
presence of a Yukawa potential (Red dots) compared with the 
modified Wannier-Stark complex spectrum (Black crosses). 
In this figure, we have chosen U = 2Er, qy ~ W^^ et Ay = 
1 /im. 



V. DISCUSSION 

As we have seen in the previous sections, the pres- 
ence of a material surface as well as the assumptions of a 
deviation from the short-range gravitational law do not 
considerably limit the lifetimes of the metastable states 
of the trap. As a consequence, the assumption consist- 
ing in considering these states as pseudo-eigenstates re- 
mains valid under these assumptions. However, up to 
now, we have not taken into account the effect of the 
Casimir-Polder interaction between the quantized elec- 
tromagnetic field in the presence of the surface and the 
atom. As we have seen in |19j, this effect is dominant 
at short distances inducing an important correction to 
the real energy levels of the atomic states. So we can 
suspect that it could have a non-negligible effect on the 
imaginary part of the spectrum as well. 

As shown in |19j . the Casimir-Polder atom-surface in- 
teraction can be computed using second-order perturba- 
tion theory on the Hamiltonian interaction term. Using 
the same notation as ^19J , the Hamiltonian describing our 
system can be expressed as 

H = Hq + Hint = H{ + iJat + H\fi/s + Hint 



i/at = hujQ\e){e\ 

2 
HwS = :: m 



(15) 






a9z+ — (l-cos(2fc/z)) 



The complete Hamiltonian is written as a sum of a term 
Hq describing the free evolution of the atomic and field 
degrees of freedom. In particular, H{ is the Hamiltonian 
of the quantum electromagnetic field, described by a set 
of modes (p, k, fc^): here p is the polarization index, tak- 
ing the values p = 1, 2 corresponding to TE and TM po- 



Far Regime 


Ay {^J.m) 


Uy 


Im(Ci') 
(units of 10'^° Er) 


rW (10^0 s) 


0.70000 


1.95007[12] 


-6.25618 


6.13778 


0.80835 


2.08378[11] 


-8.09505 


7.94185 


0.93347 


2.87073[10] 


-8.44136 


8.28161 


1.07795 


4.92177(9] 


-8.51022 


8.34916 


1.24480 


1.02282(9] 


-8.52701 


8.36564 


1.43747 


2.51801(8] 


-8.53137 


8.36991 


1.65996 


7.15759(7] 


-8.53326 


8.37177 


1.91689 


2.30712(7] 


-8.53361 


8.37211 


2.21359 


8.33762(6] 


-8.53419 


8.37268 


2.55622 


3.31836(6] 


-8.53403 


8.37252 


Near Regime 


Ay {^J.m) 


Uy 


Im(<') 
(units of 10-^°£;^) 


rW (10^° s) 


0.10000 


9.98047(11] 


-8.55444 


8.39255 


0.11659 


9.98047(11] 


-8.55682 


8.39488 


0.13594 


9.98047(11] 


-8.55945 


8.39746 


0.15849 


9.98047(11] 


-8.56119 


8.39917 


0.18479 


9.98047(11] 


-8.56046 


8.39845 


0.21544 


9.98047(11] 


-8.55548 


8.39357 


0.25119 


9.98047(11] 


-8.54216 


8.38050 


0.29286 


9.98047(11] 


-8.51542 


8.35427 


0.34146 


1.41907(10] 


-8.53314 


8.37165 


0.39811 


2.60450(9] 


-8.53385 


8.37235 



TABLE II: Table of the first ten imaginary parts of the com- 
plex energy in the first well as a function of the values of q:y 
and Ay calculated in [T9| . 



larization respectively, while k and k^ are the transverse 
and longitudinal components of the wavevector. iJat is 
the internal Hamiltonian of our two-lcvcl atom having 
ground state \g) and excited state (e) separated by a tran- 
sition frequency Wq- While 7?at is associated to the inter- 
nal atomic degrees of freedom, the term flws accounts 
for the external atomic dynamics. The interaction be- 
tween the atom and the quantum electromagnetic field 
is written here in the well-known multipolar coupling in 
dipole approximation [55], where ^i — qp {q being the 
electron's charge and p the internal atomic coordinate) 
is the quantum operator associated to the atomic electric 
dipole moment and the electric field is calculated in the 
atomic position r. 

In order to compute the correction to the complex en- 
ergy spectrum due to the Casimir-Polder effect, we have 
used a non-hermitian equivalent of the complex perturba- 
tion theory based on a redefinition of the scalar product, 
which will be called c-product, defined by [39] 



i^W 



+ 00 



dx(p{x)ip{x). 



(16) 



Using this definition, we can define the correction term 



by the same way as in standard quantum mechanics and 
obtain the expression of the second order energy correc- 
tion 



AK 



(2) 



[(v^rrvr)]' 



k^n 



E.^ 



(0) 



E 



(0) 



for a complex-scaled Hamiltonian of the form 



H" = m 



w 



(17) 



(18) 



where H^ is the non-hermitian unperturbed Hamiltonian 
and W^ is a small perturbation. We recall here that the 
first-order correction in the case of the dipolar approx- 
imation for the atom-field interaction is null due to the 
nature of the atomic dipolar operator /x. After having 



scaled the Hamiltonian (151, we are now prepared to cal- 



culate the complex correction knowing that the unper- 
turbed states can be written as 



|VS) = |5)|0.(k,fc.))Kl), 



(19) 



where n stands for the number of the well where the atom 
is initially placed (in the first Bloch band). We have to 
stress here that this state belongs to a product of three 
spaces and that the transformation Q acts only on the 
external atomic states. As a consequence, we will use the 
c-product only on the external atomic states whereas we 
keep the standard scalar product acting on the internal 
atomic states and field states. These assumptions lead 
to an expression for the second-order correction of the 
modified Wannier-Stark energy levels 



AE. 



(2) 



E 



4-CX3 



dk:,_ 



^2^yy[(^Z\HLt\e)\lp{Kh))\n,b)y 



„=i,=i El^l-El^^,+h{u; + u;o) 



(20) 



r 



where 



Hi,^^ti-I^ire^'). 



(21) 



Equation ( 20 1 can be simplified considering that in a shal- 
low trap, as in the case of FORCA-G, the Bloch bands 
above the second one are not trapped in the lattice so 
that the wavefunctions associated with these bands are 
very delocalized and the overlap with the states of the 
first band is very weak. This observation allows to re- 
strict the sum over b to its first two terms. Moreover, 



due to the extension of the modified Wannier-Stark wave- 
functions in the first band, it is sufficient to take into ac- 
count only 20 terms in the sum over n. Finally, the weak- 
ness of the imaginary part of the modified Wannier-Stark 
states as well as those of the difference of the real part 
of the Wannier-Stark spectrum with respect to the term 
h{uj+u}Q) allows us to ignore the contribution E^^ b~-^m i 



in the denominator in Eq. ( 20 ) . All these simplifications 



lead to a new expression of the correction 



A<\=/ dz 

Jo Jo 



m+lO 

d^' E [-E 

n—m— 10 p 



+00 



dk^ / d^k 



(m,l|^)(z|n,l)(m,l|z')(z>,l) 
h{uj + Wo) 



(m, l|z)(z|n, 2)(m, l|z')(z>, 2) 



H{u} + Wo) 
/ dz dz'Y^ dk, d^k V 

"'O Jo „ Jo J „-,„_in'- 



<)(k,fc„r)AWt(k,fc„r') 

m+lO 



WS^,l{z)WS,nA^') 

A^^\k,k,,r)Al'^\k,k,,r'){WSnA^)WSnA^') + WSna{z)WSnA^')) 



7i—7n — 10 



(22) 



where WSn.m{x) state for the modified Wannier-Stark and 

states taking into account the presence of the surface, ,„. 

AW(k,fc„r) = -(Op(k,fc,)|M,g •£;(re»«)|lp(k,A;,)) 



huj 

IT V 47re(] 



e^'"/'^'Ves-fp(k,fc„ze'«), 
(23) 



9 



where fp(k, fc^jz) are the mode functions of the electric 
field in the presence of a perfectly conducting mirror in 
z = gD] 

fi(k, fc^, z) = k X z sin(fc2z) 

f2(k, kz, z) = k sin(fc2z) — z — cos{kzz), 

being k = k/fc and z = (0,0,1). Unfortunately, the 



terms Jvp (k, fc 



iW^ 



iWtr 



r) and Ac^ (k, fc^jr') lead to a diver- 
gent integral when the coordinate z is complex-scaled. 
Indeed, using the definitions of the electric field given in 
[TO] , we obtain the following divergent expression for the 
correction 



m+lO 



{■■■ ' 



47r3eo 

■T^™,iWM^™,i(^') 






n— T71— 10 



(25) 



(l + ( — ) ^) sin(fc^ze'*) sin(fc^z'e*^) + [— ] cos(fc^ze*^) cos(fc;,z'e"')] } 



This divergence is fundamental because the terms in cosh 
and sinh (coming from the sin and cos functions with 
complex arguments) produce a divergence which is not 
sensitive to any ordinary regularization technique. This 
could be due to the fact that the atom is treated all along 
the calculation as a pointlike particle whereas its finite 
size should be taken into account. However, this effect is 
difficult to characterize because the scaling of the Hamil- 
tonian makes the coordinate representing the atom-wall 
distance complex so that the atomic size should be de- 
fined as a complex quantity, which is far from natural. 
So, for the time being, the calculation of the effect of the 
Casimir-Polder interaction on the lifetimes of the atomic 
states remains an open problem. 



VI. CONCLUSIONS 

In this paper we have studied the modifications to the 
lifetime of atoms trapped in an optical lattice in prox- 
imity of a surface. We have shown that the boundary 
condition introduced by the presence of the surface only 
marginally modifies the value of the ordinary Wannier- 
Stark lifetimes, leaving them almost infinite with respect 
to the duration of a measurement in a typical atomic- 
interferometry experiment such as the recently proposed 
FORCA-G. The same holds for the presence of a hypo- 
thetical Yukawa deviation from Newton's gravitational 
law. In our analysis, we have modeled the surface both 
as an infinite and a finite potential barrier, and consid- 
ered both an atom above or below the surface. 

As a natural development of our work, it could be 
interesting to investigate the behavior of the states in 
front of a more realistic surface, abandoning the assump- 



tions of perfect conduction and infinite extension. This 
is the subject of ongoing work and will be part of an 
upcoming publication. As discussed in the last section, 
another problem which remains open is the precise cal- 
culation of the perturbation of the lifetimes due to the 
Casimir-Polder effect. However, we expect that this will 
not significantly reduce the lifetime. On the contrary, 
the Casimir-Polder force is attractive, so for atoms be- 
low the surface it will tend to counteract Landau-Zener 
tunnelling thus increasing the lifetime, rather than de- 
creasing it. This assumption is made knowing that the 
Casimir-Polder potential is attractive toward the surface 
whereas the Landau-Zener effect tends to drag the atom 
away from the surface. Thus, at the moment we can as- 
sume that the Casimir-Polder effect should act on the 
lifetimes in an opposite way with respect to the Landau- 
Zener effect so that it should increase the lifetime of the 
atomic states. Finally, we point out that in this work we 
have not taken into account chemical processes between 
the atoms and the surface (atoms "sticking" to the sur- 
face) as we have considered the surface as a simple po- 
tential barrier. A more realistic analysis should include 
an additional potential describing that interaction, but 
is beyond the scope of the present paper. 
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